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The superconductivity in high temperature supercon- 
ductors ordinarily arises when doped with hetero-valent 
ions that introduce charge carriers ifiUi]. However, 
in ferropnictides, "iso-valent" doping, which is gener- 
ally beUeved not to introduce charge carriers, can in- 
duce superconductivity as well l^-ill]. Moreover, unUke 
other ferropnictides the superconducting gap in 

BaFe2(Asi vPv)2 has been found to contain nodal lines 1 14|- 
[16]. The exact nature of the "iso-valent" doping and nodal 
gap here are key open issues in building a comprehen- 
sive picture of the iron-based high temperature supercon- 
ductors 117-20]. With angle-resolved photoemission spec- 
troscopy (ARPES), we found that the phosphor substitu- 
tion in BaFe2(Asi vPv)2 induces sizable amount of holes 
into the hole Fermi surfaces, while the t/n -originated band 
is relatively intact. This overturns the previous common 
belief of "iso-valent" doping, explains why the phase di- 
agram of BaFe2(Asi vP v)2 is similar to those of the hole- 
doped compounds, and rules out theories that explain the 
nodal gap based on vanishing t/vv hole pocket. 

BaFe2(Asi_fPj^)2 is a rather unique ferropnictide as its su- 
perconductiyity is introduced by the iso-valent doping of P 
for As isi |6t]. UnHke the hetero-valent doping that alters 
the carrier concentration in Bai_ vKvFe2As2, BaFe2-.vCoi As2, 
or LaOi-^FvFeAs |2-4], the iso-valent doping is often con- 
sidered not to alter the occupation of the Fe 3d bands, as 
illustrated by the density functional theory calculations of 
BaFe2As2 and BaFe2P2 as well [6, 7]. Yet, surprisingly, it 
has a similar phase diagram just like the hetero-valent doped 
cases: with P doping, spin density wave (SDW) is suppressed 
and superconductivity (SC) emerges [6]. 

Since P anion is smaller than As anion, and thus introduces 
internal strain or distortion, i.e. chemical pressure, the su- 
perconductivity introduced by iso-valent doping is associated 
with the unprecedented pressure dependence of the supercon- 
ducting transition temperature {Tc) generally observed in iron- 
based superconductors [21-24]. In fact, it is the largest among 
all superconductors in both relative and absolute scales. For 
example, a Tc dependency of 2-4K/GPa and sometimes even 
lOK/GPa is observed in BaFe2(Asi_j;P^)2, LaOi-iF^FeAs, 
etc. [I2l[ I22II: and an increase of Tr from to above 30 K 



is observed in BaFe2As2 and FeSe under pressure | 23l 12411 . 
However, these remarkable pressure effects are still far from 
understood. Theoretically, P doping is predicted to alter the 
band structure and Fermi surface topology dramatically, con- 
sidering it changes the electron hopping terms [17,25]. Partic- 
ularly, it is predicted that the c/^^-based band would go above 
the Fermi energy {Ef), while the t/vj-based band would move 



down below Ef with P doping. Several theories further claim 
that nodes will appear in the superconducting gap when the 
d„ hole Fermi pocket disappears [17-19]. 

Figure[T]examines the dependence of the Fermi surfaces on 
the P concentration in a series of BaFe2(Asi_^Pf)2, where the 
photoemission intensity maps near Ef are shown for two ^.'s. 
The features at the zone center (F and Z) are hole pockets, 
and those at the zone corner (M and A) are electron pockets 
ll26i I281I . As P doping increases, the size of the hole pockets 
increase significantly, while the electron pockets show neg- 
ligible doping dependence. This indicates that the P doping 
could induce extra holes into the system, contradicting to the 
ordinary picture of iso-valent doping. 

To understand such extraordinary P doping effect, more de- 
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FIG. 1 : The doping dependence of the photoemission intensity 

map in BaFe2(Asi_,Pv)2. a, Photoemission intensity maps around 
r and M {k. = 0) in BaFe2(Asi_,P,)2 for x = 0, 0.15, 0.25, 0.30 re- 
spectively, named as PO, P15, P25, and P30 hereafter b, is the same 
as a, but taken around Z and A {k^ — 2n/c). Data around T, M, Z, 
and A were taken at 22 eV, 26 eV, 33 eV, and 17 eV respectively. All 
data were measured in the paramagnetic state at 150 K, 1 10 K, 75 K, 
and 40 K for PO, P15, P25 and P30 respectively to avoid the compli- 
cations from the electronic structure reconstruction in the SDW and 
superconducting states [|2^l27ll . 
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FIG. 2: The doping and polarization dependence of the photoemlssion data, a, Left side: the second derivative with respect to energy of 
the photoemission intensity taken along F-M direction with 118 eV photons for PO and P15, and 120 eV photons for P25 and P30 in the p 
polarization; Right side: photoemission intensity highlighting the fast dispersing weak features, b, is the same as the left side of a, but taken in 
the .V polarization. e,f, the corresponding momentum distribution curves (MDC's) near Ef of a,b. Note that the right part of e is in an expanded 
intensity scale. The peaks marked by black solid squares in e are contributed by the residual spectral weight of the band below Ef with the d.i 
orbital (see supplementary information for details), c, d, g, and h, are the same as a, b, e, and f respectively, but taken along Z-A direction 
with 100 eV photons. Solid lines are band dispersion, as determined by MDC peaks or the minima of the second derivative of photoemission 
intensity with respect to energy. Only half of data are overlayed with determined dispersions for a better view of the data. 



tailed band structures near Ef are examined with both s and 
p polarizations of the incoming photons. Since the low-lying 
electronic structures are mostly made of Fe-3d orbitals with 
specific even or odd spatial geometry, the photoemission in- 
tensity of the even (or odd) component of a band is only de- 
tectable with the p or i polarized light (l^ (see supplementary 
information for details). Fig.|2^ shows the photoemission data 
taken around F with the p polarization. One could resolve a 
fast dispersing a band. According to previous studies, the a 
band corresponds to the even dxz orbital [28 J. In the s polar- 
ization (Fig. I2J5), two bands, /3 and 7, could be observed. The 
/? band could be assigned to the odd dy^ orbital. Although /3 
is almost degenerate with a near Ef, it bends back towards 
Ef at about 120 meV, instead of further dispersing to higher 



binding energies. The y band is made of the d^- orbital, which 
shows much weaker intensity than /3 due to its small matrix 
element near the F point \28\. Around the Z point, however, 
one could resolve two bands in the p polarization (Fig. |2};). 
The inner band is contributed by the even component of the 
/3 band, since it shows the identical Fermi crossing and band 
dispersion with the odd component of the /3 band observed 
in the s polarization (Fig. |2}i). The outer one is the a band, 
whose Fermi momentum is enlarged significantly from F to 
Z, indicating its remarkably strong k, dispersion. On the other 
hand, the Fermi momentum of j8 moves outward only slightly 
and becomes degenerate with the y band around Z. 

Two electron-like bands, 6 and rj, are clearly distinguished 
in the s and p polarization separately around the zone cor- 
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FIG. 3: The experimental Fermi surface cross-sections in Z-F-M-A plane, a-d, Fermi surfaces of PO, PI 5, P25, and P30 respectively 
determined based on the experimental Fermi crossings (dots), e, Fermi surface of Bao.6Ko.4Fe2As2 1.27.1 . 



ner, consistent with previous studies [28]. Because of the 
weak intensity, the rj band is presented with the photoemis- 
sion intensity data in Fig. |2h,c, and the Fermi crossings of 
these two bands together with others could be traced from the 
corresponding momentum distribution curves (MDC's) near 
Ef in Fig. |2^-h. It is clear that the Fermi crossings of 6 and 
T] show negligible doping dependence, while the Fermi mo- 
menta of a, p and y increase with P doping. Particularly, the 
a Fermi crossings expand much more exceptionally than the 
others around Z. 



The different behaviors of the band structure between F and 
Z suggest a strong three-dimensional (3D) character of the 
electronic structure. Through photon energy dependence of 
the MDC's in both the s and p polarizations (see supplemen- 
tary information for details), we could determine the Fermi 
crossings of individual bands. The resulting Fermi surface 
cross-sections in the kx-k, plane are shown in Fig. |3t-d. In 
all cases, the two electron pockets 6 and rj at the zone corner 
show a weak 3D character and negligible doping dependence 
The a and p hole pockets increase from F to Z, while 
the y pockets show almost no dependence, which is consis- 
tent with the two-dimensional character of the c/ i,, orbital 12811 . 
With P doping, fi and y Fermi surfaces expand slightly, keep- 
ing their shapes basically unchanged. However, the curva- 
ture of the a Fermi surface changes significantly, which sug- 



gests a strong enhancement of the 3D character of the a band 
due to P doping. While previous quantum oscillation ^ and 
ARPES measurements [29J were focused on limited regions 
of the Brillouin zone, our data over the entire Brillouin zone 
give direct evidence for the phosphor induced hole doping. 
Moreover, the three hole pockets observed here complement 
the previous ARPES study [!()], where only two hole pockets 
were observed due to its particular experimental setup. Fur- 
thermore, our results provide compelling evidence that the d^y 
band still forms the y Fermi surface near Ef, which is incon- 
sistent with the prediction of band calculations lll7l4l9ll . Con- 
sequently, our data rule out all the theories that explain nodal 
gap in BaFe2(Asi_ ,Pj^^)2 based on the vanishing dxy Fermi sur- 
face CMl- 



The hole doped by phosphor in BaFe2(Asi_ vPj)2 is unex- 
pected. As a comparison. Fig. [3^ shows the Fermi surface 
cross-section of the optimally hole-doped Bao.6Ko.4Fe2As2 
with a Tc as high as 38 K |27], where the doped potassium 
is outside the FeAs layer and should not induce much chem- 
ical pressure in the FeAs layer Compared with the undoped 
case (Fig. [3^), the hole doping in Bao.6Ko.4Fe2As2 enlarges 
the hole pockets, and shrinks the electron pockets as expected 
in a rigid band picture, where the 6 electron pocket even di- 
minishes. On the other hand, the different behaviors of elec- 
tron and hole pockets in BaFe2(Asi_fP^)2 manifest its strong 
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FIG. 4: Summary of the band-specific carrier density, and pliase diagram of BaFcaCAsi ,Pv)2- a, The carrier density contributed by 
individual bands, and b, the total caiTier density in BaFe2(Asi_vPv)2 and Bao.6Ko.4Fe2As2. c, The phase diagram of BaFe2(Asi_,Pj-)2 with 
respect to the total carrier density. The error bars of caiTier density per Fe in a, b, and c are due to the uncertainty in estimating the 3D Fermi 
surface size. The error bars of the transition temperatures in c are obtained from the transport measurement. 



non-rigid band nature, which is most likely caused by chemi- 
cal pressure. More specifically, the y pocket is much smaller, 
and the electron pockets {?] and S) are lager in P30 than that in 
Bao fjKo 4Fe2As2. Another remarkable diff'erence between P30 
and Bao 6Ko4Fe2As2 is the a Fermi surface, which exhibits a 
strong k- dependence, and becomes the outmost hole pocket 
around Z in P30. Recent theoretical study has proposed that, 
with the increasing of the phosphor doping, the weight of d^i 
orbital could shift up towards Ef, and strongly mixed into the 
bands around Z, which could result in the nodal supercon- 
ducting gap i20ll . Since the mixture of d^i might enhance k~ 
dependence, our results seem to favor this scenario. Thus, it 
would be elucidative to study the d^i component in a, and look 
for nodes in the superconducting gap on the a Fermi surface 
around Z. We leave this for further studies. 

The charge carrier doping can be estimated from the Fermi 
surface volume. Fig. summarizes the carriers per Fe con- 
tributed by each band. The hole carriers in the a, yS, and y 
bands increase simultaneously with P doping, with the biggest 
contribution from the a band. The electron carriers in the 6 
band remain the same, while those of the t] band slightly in- 
crease with P doping. Fig. ^ plots the total net carrier per 
Fe obtained by summing up the contributions from all these 
five bands. The doping of holes is nearly proportional to the P 
doping, and one phosphor substation could induce about 0.3 
holes per iron. One possible explanation for such a hole dop- 
ing comes from the balance between bond length and anion 
electron negativity. When the bond length is shorter, an anion 
would take less electrons away from the cation of a covalent 



bond. For BaFe2As2 and BaFe2P2, the Fe-As bond length is 
2.3980 A, and Fe-P bond length is 2.2614 A respectively fl. 
The electronegativity of P (2.19 in Pauling's scale) is larger 
than that of As (2. 18) |30]. The smaller Fe-P bond length and 
the larger P electronegativity make the charge distribution on 
the Fe-P covalent bond similar to that of Fe-As. Therefore, 
holes and electrons in both systems are balanced. However, 
for light P-doped case, the P substitution could induce local 
lattice distortion and the Fe-P bond length is larger than that 
in BaFe2P2 (2.2614 A). As a result, P would take more elec- 
trons away from Fe, and the localization of electrons at phos- 
phor sites induces the hole carrier. On the other hand, when 
the doping is approaching the BaFe2P2 end, the hole doping 
would predictably decrease. 

For comparison, the carrier density distribution in 
Bao.ftKo 4Fe2As2 is shown by the pentacles in Fig.|4^,b, where 
the y band contributes the most holes into the system rather 
than the a band in BaFe2(Asi_tP.v)2- The total carrier den- 
sity in Bao 6Ko.4Fe2 As2 is 0.208 ±0.0 10 holes per Fe, as ex- 
pected from the K concentration. That is, the hole dop- 
ing in Bao.6Ko.4Fe2As2 is due to the charge transfer between 
the FeAs layer and the potassium layer, instead of the lo- 
calization of electrons in BaFe2(Asi_iPv)2- Fig. |4j; gives 
the phase diagram with respect to the total carrier per Fe in 
BaFe2(Asi_^Pf)2. It greatly resembles the K-doped case, ex- 
cept that Tc reaches its maximum at 0.09 holes per Fe, instead 
of 0.2 holes per Fe in Bai_ vKrFe2As2. It suggests that the rise 
of superconductivity in BaFe2(Asi_iPj^)2 is not a pure dop- 
ing eff'ect. Both the doping eff'ect and the chemical-pressure- 
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induced non-rigid band behavior described above should be 
considered to explain the phase diagram of "iso-valent" dop- 
ing in iron-based superconductors. 

Recently, it has been demonstrated that the chemical and 
physical pressure are rather equivalent in BaFe2(Asi_j:Pi)2, 
as the arbitrary combination of both phosphor doping and 
physical pressure can reproduce the same phase diagram 112 2|| . 
Based on our results, it is not unnatural to speculate that the 
drastic physical pressure effect in iron-based superconductors 
is likely due to the similar electronic structure changes. How- 
ever, as the physical pressure does not introduce carriers, both 
electron and hole pockets would vary simultaneously. 

To summarize, the iso-valent ionic picture is oversimpli- 
fied which could not describe the phosphor doping in iron- 
based superconductors. Our results highlight that both chem- 
ical pressure effect on the band structure, and a new route of 
carrier doping should play important roles in understanding 
the various distinct properties of BaFe2(Asi_^Pf)2- Particu- 
larly, the d„. band does not sink below Ef even at the opti- 
mally doped sample, which disproves the theories that tried 
to explain the nodal gap in BaFe2(Asi_iPv)2 based on the ab- 
sence of the dxy Fermi surface. More importantly, the phe- [18 
nomenologies in BaFe2(Asi_iPf)2 are united with those by 
carrier doping to some extent; and a possible explanation of 
the remarkable pressure effects on in iron-based supercon- 
ductors is given. 
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Methods: High quality BaFe2(Asi_^P.t)2 (x = 0, 0. 15, 0.25, 
0.30) single crystals were synthesized without flux. Ba, FeAs 
and FeP were mixed with the nominal compositions, loaded 
into an alumina tube, and then sealed into a stainless steel 
crucible under the Ar atmosphere. The entire assembly was 
heated to 1673 K and kept for 12 h or longer, and then slowly 
cooled down to 11 73 K at the rate of 4 K/h before shutting off 
the power. Shiny platelet crystals as large as 2x 2x 0.05 mm^ 
were obtained with residual resistivity ratio of about 10(see 
supplementary information for details). The P concentrations, 
X, were confirmed by an energy dispersive X-ray (EDX) anal- 
ysis. The ratio of Ba, Fe, and (AsP) is about 1:2.1:1.9 in all 
samples. The spin density wave (SDW) transition tempera- 
tures were 137 K, 100 K, and 65 K in x = (PO), x = 0.15 
(P15), and x = 0.25 (P25) samples respectively. Supercon- 
ductivity were observed in P25 and x =0.30 (P30) samples 
with = 9 K and 30 K respectively. Mixed polarization data 
were taken at the Beamline 5-4 of Stanford Synchrotron Radi- 



ation Lightsource (SSRL). Polarization-dependence data were 
taken at the SIS beamline of the Swiss Light Source (SLS). All 
the data were taken with Scienta electron analyzers, the over- 
all energy resolution was 15-20 meV at SLS or 7-10 meV at 
SSRL depending on the photon energy, and the angular reso- 
lution was 0.3 degree. The samples were cleaved in situ, and 
measured under ultra-high- vacuum of 5 x 10"' ^torr. 
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I. SAMPLE DESCRIPTION 



High quality BaFe2(Asi_vP.v)2 (x = 0, 0.15, 0.25, 0.30) single crystals were synthesized with no 
flux method. Shining platelet crystals as large as 2x 2x 0.05 mm^ were obtained as shown in the 
inset of Fig. Sla. The resistivity of BaFe2(Aso.7Po.3)2 single crystal shows a sharp superconducting 
transition at 30 K, with the transition width less than 0.5 K (Fig. Sla). The room temperature-to- 
residual resistivity ratio is « 10.4, indicating minimal disorder and impurity of the single crystal. 
The magnetic susceptibility as a function of temperature was measured in a 20 Oe magnetic field 
(Fig. Sib). A sharp drop in the zero field cooling (ZFC) susceptibility, indicating the magnetic on- 
set of superconductivity, appears at 30K, which is consistent with the zero-resistivity temperature 
in Fig. Sla. The ZFC and field cooling (FC) curves indicate a bulk superconductivity nature and 
high quality of the crystals. The phase diagram of BaFe2(Asi_;tP;t)2 is shown in Fig. S2 




Temperature(K) Temperature(K) 

FIG. SI: Electrical resistivity and magnetic susceptibility of BaFe2(Aso.7Po.3)2- a, Temperature de- 
pendence of electrical resistivity. The inset is the photography of the BaFe2(Aso.7Po.3)2 single crystal, b, 
Temperature dependence of magnetic susceptibility in a 20 Oe magnetic field. 
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FIG. S2: Phase diagram of BaFe2(Asi-vPA^)2 with respect to the P doping. The P concentrations, x, were 
confirmed by an energy dispersive X-ray (EDX) analysis. With P doping, the spin density wave (SDW) is 
suppressed and the superconductivity (SC) emerges. The SDW transition temperatures are 137 K, 100 K, 
and 65 K in X = (PO), x = 0.15 (P15), and x = 0.25 (P25) samples respectively. Superconductivity is 
observed in P25, x =0.30 (P30), and x =0.37 samples with Tc = 9 K, 30 K, and 24.5 K respectively. 
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II. EXPERIMENTAL SETUP 



The polarization-sensitivity of the orbitals in angle-resolved photoemission spectroscopy 
(ARPES) is a powerful tool to identify the orbital characters of band structure dli |2|]. The ma- 
trix element of the photoemission process can be described by 

|M)/cxK0j|e.r|0f)p 

, where s is the unit vector of the electric field of the light . For high kinetic-energy photoelec- 
trons, the final-state wavefunction 0^. can be approximated by a plane-wave state e*''^ with k in 
the mirror plane as plotted in Fig. S3c. Consequently, it is always even with respect to the mirror 
plane. For the p (or s) experimental geometry in Fig. S3c, because e is parallel (or perpendicular) 
to the mirror plane, e ■ r is even (or odd). Therefore, to give finite photoemission matrix element, 
i.e. to be observable, the initial state has to be even (or odd) in the p (or s) geometry. 

Considering the spatial symmetry of the 3d orbitals (Fig. S3d), when the analyzer slit is along 
the high symmetry direction of the sample, the photoemission signal of certain orbitals would 
appear or disappear by specifying the polarization directions as summarized in Table |Sl] For 
example, with respect to the mirror plane formed by direction #1 and sample surface normal (or 
the xz plane), the even orbitals {d^y, d,i, and d^2_yi) and the odd orbitals (J^)' and dyy) could be 
only observed in the p and s geometry respectively. Note that, d^. and dy^ are not symmetric with 
respect to the mirror plane defined by direction #2 and surface normal, thus could be observed in 
both the p and s geometries. 



TABLE S 1 : The possibility to detect 3<i orbitals along two high symmetry directions in the p and s 
geometry by polarization-dependent APRES. 
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FIG. S3: The experimental setup and definitions, a, The illustration of the Brillouin zone of BaFe2As2. b, 
Two-dimensional plot of the simplified Brillouin zone (solid line) and the unfolded Brillouin zone (dashed 
line), c, Experimental setup for polarization-dependent ARPES. For the p (or s) experimental geometry, 
the electric field direction of the incident photons Sp (or s^) is parallel (or perpendicular) to the mirror plane 
defined by the analyzer slit and the sample surface normal, d, Illustration of the spatial symmetry of the 3d 
orbitals with respect to the mirror plane formed by surface normal and cut #1 in panel b, i.e. the xz plane, 
e, A typical orbital assignment of bands of iron pnictide as calculated in Ref. [3]. 
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III. d^2 ORBITAL 



Besides of the fast dispersing a band, another peak-like feature could be observed in the mo- 
mentum distribution curves (MDC's) in the inset of Fig. S4a, as pointed out by an arrow. Note 
that, the d^2 orbital contributes to one band (^), which is 200 meV below Ep (Fig. S3e) M\- This 
band disperses towards the a band, resulting in an enhancement of the spectral weight at the band 
crossings (Figs. S4a, b). The resulted intense spectral weight would extend some residual weight 
towards Ef (Fig. S4c) due to its broad peak width, which forms the non-dispersing peaks in MDC's 
as marked by the solid squares in Fig. S4d. Therefore, the peak-like feature (^) observed in the 
p polarization at is not the Fermi crossing of band, but the residual spectral weight of ^ well 
below Ef. 
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FIG. S4: The photoemission data of P30 near Ef taken in the p polarization, a, b, The photoemis- 
sion intensity {I{k,a))) and its second derivative with respect to energy {d^I{k,aj)/daj^) taken with 120 eV 
photons around T in the p polarization, respectively. Solid lines are band dispersion, as determined by mo- 
mentum distribution curve (MDC) peaks or the minima of the second derivative of photoemission intensity 
with respect to energy. The inset in a is the MDC at Ef. c, d. The energy distribution curves (EDC's) and 
MDC's of the data in a, respectively. 
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IV. PHOTON ENERGY DEPENDENCE 



By changing the photon energy in ARPES experiment, one can probe the different values, 
according to the following formula [|4|] : 

K = ^[(hv - O)cos20 + VoV'^ 

where m and h are the electron mass and Plancks constant, v is the photon energy, O is the work 
function, and Vq is the inner potential of the sample, which could be determined experimentally. 
Through matching the periodicity of the photon energy dependence to the high symmetry points 
of the Brillouin zone along direction, the inner potential (Vq) is set to 16 eV for calculating the 
fc^'s in our experiment. 

Fig. S5 shows the photon energy dependence of the MDC's near Ep. The Fermi crossings of 
different bands could be traced separately in the s and p polarization. The peak positions of the 
MDC's are marked by combining the highest intensity and the center of mass, which is a reliable 
way and is practiced commonly in ARPES community. Note that, while the a band move outwards 
from r to Z, the even J^z orbital in the a band could mix into the /3 band. As shown in Figs. S5e-h, 
when the L moves away from F, two peaks could be identified corresponding to the a and /3 bands 
in the p polarization. The intensity transfer between the MDC peaks of a and /3 observed in the p 
polarization is due to the k^ dependence of the even orbital mixing between the a and band. 
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FIG. S5: The photon energy and doping dependence of the MDC's near Ef in the s and p polariza- 
tions, a-d, The photon energy dependence of the MDC's near E/? taken in the s polarization in PO, P15, P25, 
and P30 samples, respectively, e-h, are the same as a-d, but taken in the p polarization. Note that the right 
parts of e-h are in an expanded intensity scale. The peaks marked by black solid squares are contributed by 
the residual spectral weight of the band below Ef with the d^i orbital 
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